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ABSTRACT 


A theoretical behavioural study of composite beams in 
negative bending is presented, with emphasis on the effect of 
longitudinal slab reinforcement. A method for determining the 
moment producing local buckling of the compression flange is 
developed and a computer program for the plotting of the moment « 
curvature and related relationships is presented. Analytical 
results are compared with the results of tests previously 
conducted at the University of Alberta. 

Two approaches are used in obtaining the flexural 
behaviour relationships; one based on an idealized stress-strain 
diagram consisting of elastic, plastic and strain-hardening 
regions, the other based on a modified stress-strain diagram 
consisting of elastic and strain-hardening regions. Restraint 
afforded the flange by the web is based on the location of the 
neutral axis. An upper limit is placed on the effectiveness of 
the longitudinal slab reinforcement. 

The analysis based on the modified stress-strain diagram, 
in which the two regions are defined by an elastic modulus, E, 
and a strain-hardening modulus, ESe> provides a satisfactory 
method for determining the moment at local buckling for 
composite beams in which strain-hardening strains are reached in 


the compression flange. Although there is a considerable 
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variation between theoretical and experimental values of curvature 
at ultimate moment, it is considered unlikely that any theoretical 


analysis will provide any better agreement. 
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CHAPTER I 
INTRODUCTION 


lel INTRODUCTORY REMARKS 

In order to realize the full advantage of composite construction 
applied to flexural members it is necessary to develop continuity. The 
behaviour of composite beams in resisting positive moment is well 
defined; however, the behaviour of composite beams under negative 
moment is not as fully understood. In negative moment regions 
longitudinal slab reinforcement may be introduced to increase the 
moment capacity since the concrete slab is ineffective in tension. Test 
results indicate that this reinforcement effectively increases the 
negative moment resistance and current design specifications permit full 
contribution of longitudinal reinforcement. 

In the design of structural steel flexural members, care is 
taken to prevent premature local buckling of the compression elements. 
Local buckling must also be considered in composite beams, particularly 
in continuous spans. Conditions in positive moment regions present 
little problem as the presence of the concrete slab in the compression 
zone tends to reduce the depth of the steel section in compression, 
thereby reducing the probability of local buckling. Furthermore, the 
slab tends to stiffen the compression flange. In negative bending, 
however, the concrete slab is connected to the tension flange of the 


steel section. If longitudinal reinforcement is present, most of the 
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steel section is in compression. Also, in negative bending, the 
compression flange is remote from the slab. Therefore, local buckling 
is more significant in negative moment regions. 

Recent experimental investigations have provided information 


concerning the buckling behaviour of composite beams in negative bending. 


Tee RESEARCH AT THE UNIVERSITY OF CAMBRIDGE 

Tests simulating the negative moment region of a continuous 
beam were carried out by van Dalen in 1966 on simple span composite 
beams. The modes of failure observed included crushing of the bottom 
surface of the concrete slab, longitudinal splitting of the slab, shear 
failure in the slab, buckling of the compression flange and buckling 
of the web of the steel section. 

Results of this investigation indicated that the resistance 
to bending in a negative moment region, with longitudinal slab 
reinforcement present, continues to increase with curvature until 
buckling occurs in the steel section. This was true provided adequate 
shear connection and transverse slab reinforcement were provided. 
Cracking of the concrete slab, due to tensile forces present, did not 
affect this behaviour. It was concluded that, using the assumption of 
complete interaction and plane sections, reliable predictions of the 
observed longitudinal strains, curvature and rotations, at all loads 


up to the ultimate load, may be obtained. 


Es RESEARCH AT THE UNIVERSITY OF ALBERTA 
Three experimental investigations have been conducted in the 


Department of Civil Engineering at the University of Alberta concerning 
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3. 
the behaviour of composite beams in negative bending. The beams tested 
consisted of rolled steel sections connected by headed stud shear 
connectors to reinforced cast-in-place concrete slabs. 

In 1968 Piepgrass (1) tested six specimens representing an 
isolated negative moment region. The main variables were slab width 
and the amount of longitudinal slab reinforcement. Four of the six 
beams were coverplated. One of the remaining beams was coverplated 
after the initial test and retested. The failure mode of the beams 
without coverplates was local buckling of the flange and web. The 
beams with coverplates failed by lateral buckling. It was concluded 
that the ultimate capacity of the composite section increased as the 
longitudinal slab reinforcement was increased. It was also indicated 
that a steel section defined a compact section under the provisions 
of current specifications had significantly reduced rotation capacity 
in a composite section as compared with that of the plane beam. This 
reduced rotation capacity might be of concern in plastic design in 
relation to the effect on the development of a mechanism. 

Davison (2) in 1969 studied the effect of longitudinal slab 
reinforcement on the rotation capacity of composite sections in 
negative bending. Seven beams with varying slab widths were tested. 
Failure of the beams occurred by local buckling of the web, local 
buckling of the compression flange or a combination of the two, 

The addition of longitudinal slab reinforcement was found 
to significantly increase the ultimate negative moment capacity. The 
increase, however, was not proportional to the increase in simple plastic 


moment values. As the longitudinal slab reinforcement increased, the 
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rotation capacity decreased. The width of the concrete slab had no 
Significant effect on the ultimate moment capacity or rotation 
capacity. The longitudinal slab reinforcement reached yield prior 
to the attainment of ultimate moment. 

In a third study by Lever (3) in 1970, nine beams were tested 
in negative bending. The main objective was to study the effect of 
varying the size of the steel section and to further study the effects 
of varying the amount of longitudinal slab reinforcement. The beams 
failed by local buckling of the compression flange and web, either 
separately or in combination. 

The major conclusions drawn from the three investigations 
carried out at the University of Alberta are as follows: 

1. The addition of longitudinal slab reinforcement can result 
in significant increases in the ultimate moment capacity 
of composite beams in negative bending. 

2. The increase in ultimate moment is not directly 
proportional to the increase in the theoretical plastic 
moment. 

on For a given steel section, a reduction in the rotation 
capacity of a negative plastic hinge occurs with an 
increase in the amount of longitudinal slab reinforcement. 

4. Slip between the concrete slab and steel section increases 
with load. The slip becomes concentrated in the region 
of the beam nearest to the local buckle after buckling 


has occurred. 
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5. Slip between the concrete slab and the steel section 
Significantly affects the strains in the longitudinal 
reinforcement. In all tests, strains at ultimate load 
reached the yield strain, but were considerably less 
than the value determined by extending the strain profiles 
from the steel section to the level of longitudinal 
reinforcement. 

6. Cracking behaviour of the concrete slab is influenced 
by the amount of longitudinal slab reinforcement, the 
Spacing of the transverse reinforcement, and the 
transverse curvature of the slab. A few wide cracks 
develop for low amounts of reinforcement. Increasing 
the amount of slab reinforcement results in an increase 


in the number of cracks and a decrease in the crack width. 


1.4 SCOPE OF THE PRESENT INVESTIGATION 

The primary objectives of the present investigation are to 
obtain theoretical values for local buckling moments and to develop 
theoretical moment-curvature relationships for composite beams in 


negative bending. 
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CHAPTER II 
LOCAL BUCKLING OF STEEL BEAMS 


‘ae PLATE BUCKLING IN THE ELASTIC RANGE 

Local buckling of the compression flange of an "I" section, 
due to axial load or bending, is essentially the buckling of a plate 
Subjected to uniformly distributed compressive stresses. The correct 
differential equation of a buckled plate was first developed by Navier 
in 1823. The solution of this differential equation is credited to 
Bryan in 1888. For a simply supported rectangular plate with edge 
compression in one direction, the differential equation has the form 


6 4 4 4 
Et ow A 20 WW 1 ) WwW 4 ei t 3 2y (221) 
fai ay \saee a XS) cee ve Ox? 


where fy is the edge compression stress as shown in FIGURE 2-1. The 
deflected shape may be assumed in the form of a series, each term of 


which satisfies the boundary conditions. 
c8 o8 


ay eZ 


sin 
Tooele Ann d 


~ singe (2-2) 


From EQUATIONS (2-1) and (2-2) the critical stress is found to be 


Deaaunae cay ale bame cats ee 
in which the quantities m and n represent the number of half-waves 


{nto which the plate buckles in the x and y directions. The term 


(nb + 2.0")? is defined as the plate buckling coefficient "k". 
a n 


Regardless of the boundary conditions, the critical stress for any 
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Te 
rectangular plate subjected to edge compression in one direction may 


be expressed as 


z Pa Sed ba ek 
Pee kat aatS (2-4) 


in which k, the plate buckling coefficient depends on the boundary 
conditions. 

For a comparatively long steel plate, with boundary conditions 
representative of the flange of an "I" section (y = 0 fixed and 
y = b free) the value of k is 1.328 and 

epee 2000 Gabe 


From which 

(D/ tcp 29190) 7Fam (2-5) 
where b equals half the total flange width. Since the buckling is 
elastic, Fo, is limited to the yield stress of the material, Fy. 
EQUATION (2-5) may then be written as 

(b/t), = 190/vF, | (2-6) 
For G40.12 steel, with a yield stress of 44ksi. the critical b/t 
ratio is 29. The most slender flanged "I" or "WF" section available 
th Cadhadas a OWF 15.5, has a b/t ratio of 22.3. This implies, there- 
fore, that all available "I" and "WF" sections will not buckle locally 


in the elastic range. 


ne POST-ELASTIC PLATE BUCKLING 


One of the earliest solutions to the problem of post-elastic 


plate buckling, by Bleich (4), considers yielding but not strain- 


hardening. The plate buckling stress is given simply and conservatively 
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where y = E,/E. The solution to EQUATION (2-7) is normally of the 
form Fp/vy being equal to a value. From tabulated values it is then 
possible to obtain Pree As the ratio F.p/vy increases, an becomes 
asymptotic to the yield stress of the material, Fy 
The material is heterogeneous during the yielding process 
with yielding taking place in localized slip bands in which the strain 
jumps from its value at the elastic limit to that at the outset of 
strain-hardening. The material is again homogeneous when it has all 
been strained to the strain-hardening range. | 
On the basis of orthotropic behaviour, Haaijer (5) derived 
expressions for the buckling strength of steel plates which predict 
the strain at buckling, for strains exceeding the elastic limit. He 
considered a rectangular plate subjected to uniform compressive stresses 
in one direction large enough to produce strain-hardening. Due to 
strain-hardening changes may occur in the values of the tangent moduli 


EY and Ee . POlsson Ss acio Hy and Uys and shear modulus G- 


Vic 
Based on relations between increments of strain and stress, 
the bending and twisting moments in terms of the deflection, w, can 


be written as 
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where I = Vl, 


The equilibrium of the bent plate can be expressed by the differential 


equation 
3° uw) 3*w SUI ese Beer Ka 
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This differential equation, (2-8), can be solved if H? is assumed equal 
to D,Dys an assumption originally made by Bleich. 

The value of fy at which bifurcation of the initially perfect- 
ly plane plate occurs can be determined from EQUATION (2-8). This 
condition can also be expressed in terms of work, with the additional 
work done by the external forces due to bending of the plate being 
equal to the change in internal energy in the plate. This results 


in the integral equation 
fie wy? Zn \? OER 
Salles! dx dy J {[o(3) (Duy + Dyx) (SI) 


Be b 
: aee( so dx dy (2-9) 


If external restraints exist along the plate edges, additional terms, 
representing the work done by the restraints, must be introduced in 


the right-hand side of EQUATION (2-9). Approximate solutions can be 
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10. 
obtained by assuming the shape of the deflected surface. 
The deflection surface for a rectangular plate loaded in 
compression in the x direction and with one unloaded edge free and 
the other restrained against rotation (conditions similar to those 


occurring in the flange of an "I" section) may be considered as 


W = E 7 of) + a, + 2) }] sing (2-10) 


The ratio 8 = B/A depends on the degree of restraint. 


Substituting (2-10) in (2-9) and integrating gives 
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For the most conservative case, that represented by B = 0 (one 
unloaded edge free and the other hinged), the equation for the 


critical stress becomes 
ea Gt Oey Alon | a 
F.. 2 (| [—<C) fei A) 
a ee ‘ 
Haaijer suggested the following properties based upon theoretical 


and experimental investigations 
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According to Haaijer, the yield stress of the material has little 
effect on the above values. 

Lay (6) proposed a solution to the problem of local buckling 
which eliminates some of the empiricism present in Haaijer's solution. 
It is assumed that the flange is restrained against buckling by a 
torsion spring (FIGURE 2-2a) representing the effect of the web. For 


this-case the torsional buckling equation is 


2 2 
peal cue fa Sigel Pigs) ) 
Ya eet L st “w s\ntt (2-12) 
where n = an integer; L/n = half wavelength of the buckle; 


i is the warping constant; Ky is the uniform torsion constant; and 


K, is the spring constant of the rotational spring. In EQUATION (2-12) 
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w  16\ 144 
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The spring constant, k,, can be assessed by considering the 
resistance of the web to a rotation of the compression flange. 
Assuming a deformed shape as shown FIGURE 2-2b with no rotation of 
the tension flange and no relative deflection of the two sides of the 
compression flange, the value of Ke is 
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where une half wavelength of local buckle. 


Substituting values for Es¢, Iw and ke results in 
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Based on EQUATION (2-14) Lay and Galambos (7) proposed a 
theory for local buckling of the compression flange of beams under 
moment gradient. From FIGURE 2-3 the length of the yielded region, 
etL, iS given by similar triangles as 
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It is necessary for tL to be of sufficient length, 22) for a Jocal 
buckling wave to form, where his given by EQUATION (2-14). Fora 
beam under moment gradient, one end of the yielded region is adjacent 
to the elastic part which is relatively stiff, while the other end 

is adjacent to a load point or connection as shown in FIGURE 2-4. 
Thus, relatively stiff end restraints are provided and a full wave- 


length is required for the local buckle to form. Therefore, 
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Generally, the ultimate moment obtained from EQUATION (2-15) will be 


conservative, 


(ages COMPACT SECTIONS 
The limitation on the proportions of the flange based on 
local buckling considerations can be derived from EQUATIONS (2-11) 


and (2-13). For a relatively long plate EQUATION (2-11) may be 


written as 
ee i / by aaiGe) 
or 
COVA re ev Gee (2-16) 


UponesubStIcCUETON Of 1 = Ul3, Ge = E/2.6, E = 29500 ksi and h = 40 
in EQUATION (2-13), G, = 2610 ksi. Equating the critical stress to 
Fy results in 

AS ek OCA (2-17) 
where b is half the total flange width. EQUATION (2-17) implies that 
TOE ab2/ vi local buckling of the flange will not occur until 


the flange is fully yielded. 
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On the basis of torsional buckling relationships, Lay (6) 


proposed that for the flange to be fully yielded before local 


buckling occurs. 


vie: 12.64 


For average values of material properties, this reduces to 


b < 54/vF, (2-18) 
4 , 


where Fy is considered to be POE oc 


CSA Standard S16-1969, Steel Structures for Buildings (8) 
requires that flanges subject to compression involving plastic hinge 
rotation under ultimate load have a width-thickness ratio no greater 


than 54/VFy. 


2.4 IN-PLANE BEHAVIOUR OF STEEL BEAMS 

FIGURES 2-5 and 2-6 illustrate two possible cases of loading, 
one resulting in a region of uniform maximum moment, the other 
producing a moment gradient throughout the span. For both types of 
loading, the moment-curvature relationship is linear up to the moment 


at first yield of the compression flange (M Above My» curvature 


0): 
increases more rapidly with increase in moment as yielding penetrates 
through the flange into the web with the moment approaching the plastic 
moment value. At this point the moment-curvature response for the 

two cases of loading begin to differ. For the beam under uniform 


moment, the moment becomes asymptotic to the plastic moment value 


Oe The compression flange is now completely yielded over an 
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appreciable length and it is possible for a local buckle to form. 
For the moment gradient case, moment values greater than 

My are attained due to the presence of strain-hardening. As the 
moment rises above Mp» the yielded zone of the compression flange 


lengthens until it reaches a length sufficient to cause local 


buckling. 


pS) STRESS-STRAIN RELATIONSHIPS 

In establishing the moment=curvature relationships up to 
ultimate moment conditions, it is necessary to employ the stress- 
Strain relationships for the material. The basic stress-strain 
diagram for structural carbon steel is shown in FIGURE 2-7. Three 
basic regions are shown up to ultimate stress: (1) elastic, (2) 
plastic and (3) strain-hardening. The strain-hardening modulus, 
E542 1s defined as shown in FIGURE 2-8. As the actual value of Ect 
varies non-linearly with increase in strain (FIGURE 2-7), use of the 
idealized value will result in calculated values of curvature being 
less than actual values if the ultimate moment is based on 
EQUATION (2-15). The effect of residual stresses does not appear 
in the stress-strain diagram. 

A cross-section containing residual stresses, has certain 
fibers yielding before others when load is applied to the section. 
Therefore the section is no longer homogeneous. This results in 
some difficulty in obtaining theoretical moment-curvature relation- 
ships. 


The effect of residual stress and strain-hardening, in a 
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Flange subjected to moment gradient, may be approximated by employing 
a modified stress-strain diagram as shown in FIGURE 2-9. The 
modified relationship attempts to describe the average behaviour of 
the section. The elastic range is followed by a range in which the 
material undergoes the effects of yielding and strain-hardening. 
The modified modulus Bet is less than the strain-hardening modulus 
Es¢ and its value may be obtained from the moment-curvature relation- 
ships established by tests. For a beam under moment gradient 
(FIGURE 2-6), the slope of the moment-curvature diagram above M 


p 
equals bees where I is the moment of inertia of the section. 
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FIGURE 2-1 


BUCKLED PLATE 
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(a) Flange Model 


d= 2t 


(b) Assumed Deformed Shape 


FIGURE 2-2 FLANGE BUCKLING MODEL 
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FIGURE 2-3 BENDING MOMENT DIAGRAM 
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FIGURE 2-4 SHAPE OF LOCAL BUCKLE 


105 


Connection or 
Load Point 


| sips . 
a 
, 


Rpt ii he ° 


12 agi Poanan: 
Surin? hes. 


LINQ JAO0U FO. ME cms galas 


20. 


Local Buckling 


Curvature 


FIGURE 2-5 BEHAVIOUR OF A LATERALLY SUPPORTED BEAM UNDER UNIFORM MOMENT 
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FIGURE 2-6 BEHAVIOUR OF A LATERALLY SUPPORTED BEAM UNDER MOMENT GRADIENT 
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FIGURE 2-7 STRESS-STRAIN DIAGRAM FOR STRUCTURAL CARBON STEEL 
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FIGURE 2-8 IDEALIZED STRESS-STRAIN DIAGRAM 
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FIGURE 2-9 MODIFIED STRESS-STRAIN DIAGRAM 
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CHAPTER III 
THEORETICAL ANALYSIS OF COMPOSITE BEAMS IN NEGATIVE BENDING 


Sa) BEHAVIOUR OF COMPOSITE BEAMS IN NEGATIVE BENDING 

FIGURE 3-1 illustrates the general behaviour in terms of the 
moment-curvature relationship for composite beams tested in negative 
bending at the University of Alberta. Initially curvature increases 
in direct proportion to the moment. Beyond a moment value, Mo? 
curvature increases more rapidly as the midspan moment approaches the 
plastic moment value, Wee Further increases in moment is produced by 
Strain-hardening. Local buckling of the web or compression flange 


at ultimate moment, M » results in unloading of the section with 


ult 
further increase in curvature. The moment-curvature relationship for 
composite beams in negative bending is, therefore, similar to that 
illustrated previously for plane steel beams under moment gradient. 
Since in negative bending, the compression flange is not in 

contact with the concrete slab, it behaves in a manner similar to the 
compression flange for a plane "I" or "WF" section. The concrete 
transfers shear from the shear connectors to the longitudinal slab 
reinforcement, but does not contribute to the flexural strength of 
the composite section except for very low loads. The location of the 
neutral axis depends on the amount of longitudinal slab reinforcement. 
As loading is increased, jnelastic behaviour is initiated in the 


compression flange and the neutral axis moves from the centroid of 
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the section towards the tension flange. 


oe BASIC ASSUMPTIONS FOR THEORETICAL ANALYSIS 

It is assumed that a linear strain distribution exists over 
the entire section as shown in FIGURE 3-2. To account for possible 
Slip between the steel section and the concrete slab, it is assumed 
that the longitudinal slab reinforcement yields but does not strain 
sufficiently to strain-harden. From experimental results (2) and 
(3) it has been established that the longitudinal slab reinforcement 
yields prior to ultimate conditions if the number of shear connectors= 
is sufficient to develop the ultimate moment. 

From experimental results (3) it appears that there is an 
upper limit to the ultimate moment regardless of the amount of 
longitudinal slab reinforcement. Additional area of longitudinal 
reinforcement beyond the area required to place the neutral axis in 
the tension flange results in no appreciable increase in ultimate 
moment. While the section is elastic the neutral axis is located at 
the centroid of the section. The initiation of inelastic action in 
the compression flange, due to moment values greater than es results 
in the neutral axis moving towards the tension flange. However, once 
strain-hardening of the compression flange begins the neutral axis 
migrates towards the compression flange. Therefore, the neutral axis 
is located closest to the tension flange for a moment value between 


M andM.. Assuming that the critical moment value is ee the amount 


of slab reinforcement which effectively contributes to the flexural 


strength of the composite section is, from FIGURE 3-3 
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(hw) (F) 
(ea) 
where Fw = yield stress of web, Fs = yield stress of longitudinal 


Slab reinforcement, hw = area of web and A. = area of longitudinal 
Slab reinforcement. In the theoretical analysis, if the neutral axis 
at My falls in the tension flange, the amount of longitudinal slab 
reinforcement is reduced to the amount corresponding to a neutral 
axis location at the inner face of the tension flange. 

The expression for the web spring constant, k., given by 
Lay (6) may be modified to account for the strain distribution present 
in the web. Lay assumed that the web had fully yielded under the 
longitudinal stresses present. When local buckling of the compression 
flange occurs, only the compressed part of the web is affected as 
Shown in FIGURE 3-4. As the tension zone of the web is stiff relative 
to the compression zone of the web, it is reasonable to assume that 
the compressed portion of the web acts as a propped cantilever as shown 
in FIGURE 3-4. Assuming the plate to be fixed at the location of the 


neutral axis, the spring constant expressing the restraint offered 


by the web is 


k. = (3-2) 


where z is the distance to the neutral axis from the compression 


flange, w is the thickness of the web and G; is the inelastic shear 


modulus. Thus, k. is actually a function of load for composite 


S 


sections. However, for simplicity, z is assumed to correspond 
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to the neutral axis location at moment Me 


es PROCEDURE FOR OBTAINING THEORETICAL MOMENT-CURVATURE 
RELATIONSHIPS 
The moment at first yield of the compression flange, Mo: is 
required when solving for Mizz by using EQUATION (2-15) 
M 
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The value of Mj can be obtained from equilibrium considerations, with 
the amount of longitudinal slab reinforcement equal to or less than 
the limiting value given by EQUATION (3-1), and the compression flange 
in the yielded state. An alternate method for determining Maat 
involves the calculation of Hs and the length of compression flange 
involved in the formation of a local buckie. FIGURE 3-5 shows how 

by similar triangles it is then possible to calculate the ultimate 
moment. 

In order to establish one value for the moment-curvature 
relationship, values for curvature and neutral axis location are 
assumed. It is then possible to determine the strain distribution 
over the cross-section. From the stress-strain relation the stress 
distribution and the resulting forces acting on the section are 
determined. Summing the forces will indicate whether or not the 
section is in equilibrium. If the section is not in equilibrium it 
is necessary to relocate the neutral axis and recalculate the strain 


distribution, related stresses and forces until equilibrium of the 
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section is obtained. Once the section is in equilibrium it is a 
simple matter to obtain the resisting moment corresponding to the 
assumed curvature. The complete moment-curvature relationship is 
obtained by repeating the above procedure until the calculated moment 
equals the moment at which local buckling of the compression flange 


occurs. 


3.4 TOTAL ROTATION AND CENTRELINE DEFLECTION 
Once the moment-curvature relationship is established, the 
total rotation and centreline deflection of the beam may be obtained 
by the moment-area method. Two ranges must be considered: (1) the 
entire beam in the elastic range, (2) a portion of the beam in the 
post-elastic range. FIGURE 3-6 shows the bending moment and curvature 
diagrams for the two conditions for the type of loading used in the 
University of Alberta tests. 
A) Entire span in elastic range ($ < $5) 
1) Total rotation: 
or 
2) Centreline deflection: 
aa Ged (sik 407) 
B) Post-elastic case, idealized stress-strain diagram ($ > $,) 


1) Total rotation 
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C) Post-elastic case, modified stress-strain diagram (o > bo) 
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S09 COMPUTER PROGRAM 


A computer program was developed for the calculation of the 
ultimate moment based on local buckling and woiteseece moment-curvature, 
moment-rotation and load-deflection relationships. The program, which 
is described in Appendix A, also computes the simple plastic moment 


and the limiting area of longitudinal slab reinforcement. 
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(c) TYPICAL MOMENT CURVATURE DIAGRAM 


FIGURE 3-1 UNIVERSITY OF ALBERTA TESTS OF COMPOSITE BEAMS IN 
NEGATIVE BENDING 
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FIGURE 3-2 CROSS-SECTION OF COMPOSITE BEAM IN NEGATIVE BENDING 
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FIGURE 3-3 STRESS DISTRIBUTION CORRESPONDING TO PLASTIC MOMENT WHEN 
NEUTRAL AXIS IS AT WEB-FLANGE JUNCTION 
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FIGURE 3-4 ASSUMED DEFORMED SHAPE OF CROSS-SECTION FOR 
DETERMINING SPRING CONSTANT, kK. 
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CHAPTER IV 
RESULTS 


4.1 INTRODUCTION 

Experimental results of tests conducted on composite beams 
in negative bending at the University of Alberta have been presented 
by previous investigators (2), (3). For comparison purposes, exper- 
imental moment-curvature relationships are presented... Theoretical 


results are presented in graphical and tabular form. 


4.2 MATERIAL PROPERTIES 

TABLE 4.6 represents results of tensile tests performed on 
web and flange coupons by previous investigators (2), (3). The yield 
stress and ultimate stress for the reinforcing bars used in the test 


beams is also shown. 


53 GENERAL BEHAVIOUR 


The experimental moment-curvature relationships are shown 
in FIGURES 4-1 and 4-2, while the theoretical relationships, as plotted 
on the Calcomp plotter, appear in FIGURES 4-3 to 4-6. The theoretical 
relationship between midspan bending moment and the total rotation of 
the beam is shown in FIGURES 4-7 and 4-8. Load-deflection relation- 
ships are shown in FIGURES 4-9 and 4-10. The location of the neutral 
axis, with respect to the centroid of the steel section, is plotted 


against midspan bending moment in FIGURES 4-11 and 4-12. 
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4.4 ULTIMATE LOAD CONDITIONS 

TABLE 4.1 is a summary of experimental results for tests 
conducted at the University of Alberta by Davison (2) and Lever (1). 
TABLES 4.2, 4.3, 4.4 and 4.5 summarize the theoretical study. 
FIGURE 4-13, derived from test data, is a plot of the ratio of ultimate 
moment of the composite section to ultimate moment of the steel 
section alone (Myj+/Mgi¢) versus the total tensile force in the 
longitudinal slab reinforcement at yield (Ac Fys). Corresponding 
values derived from the theoretical study are shown in FIGURE 4-14. 
FIGURE 4-15 indicates the variation in the ratio of the simple plastic 
moment for the composite section to the simple plastic moment for 
the corresponding plain steel section (Mp/Mp) with total tensile force 
in the longitudinal slab reinforcement (As PIE) The effect of 
varying the b/t ratio on the ultimate capacity of composite beams with 


constant A,/Ayp ratios is shown in FIGURE 4-16. 
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Average Slab Reinforcement Data 


mere eee ere rn CE EE 


size Yield Stress Ultimate Stress 
(ksi) (ksi) 

#3 52.4 70.02 

#4 uae 78.6 

#9 3) bal 79.4 


TABLE 4.6 ACTUAL MATERIAL PROPERTIES 


LP 


Saher S5kSeget —? 2eR8 


hy ed) 


Petty 


naa 

ai10. 3 
Par ¢ P 
OPI is 
pean oO 
eto |G 


DOLE 
ASS 


VES 


; LA 
piigess 
re 


—_ jrigasrwp 
see ager 


at 


i, de Le, pw, 


a’ 


= 


moe + bp 
ee eet 
cue) vy tae (en 
i > 


<a PAS ei 


ath 408 


Ar ee 
= 7 
_ 6 
7 


40. 


WYYOVIG JYNLVAUNI-LNSAWOW IWLNSWIYddxX4 I-t JYndI4 


( 2-OTX) HONI/SNYUTOBY-SYNLYAYNI 


Oh*T Gent UG. 08°0 03d Oh*O Uea0 00°O 


SesMica 


8I Wd 


W 
ee ZI Wad 
€l wvag pl WW3d 


00°081 00°021 00°09 00° 
SdIM’ LI-LNSAWOW 


00 °Ohe 


00 *O0€ 


r™ivn 


e 
‘ 


fer 


= 
per 
4) 
fons 
£3 
is 
3! 
wal 
< 


41. 


WYYODVIC SYUNLVAYNI-LNSAWOW TWINSWIYSdX4d ?2-b SJYNOIA 


( 2-O1X) HONT/SNYTObY-SYENLBAYND 


1270 81°O elsg 280 60°O 90°0 £0°0O 
S*OLGel 
lp WWa9 
2b WHI 
ep wag 


00°O 


00°SL 00°0S 00°S2 00°0 
SdIM° L3-LNSAWOW 


00°O01 


G0 “Se 


aa ave 0 


a 


-~ * 


42. 


(WYYOVIG NIVYLS-SS3YLS G3ZI1W3QI) 
SATYAS 9EAMCT YOS WVYOVIC JYNLVAYNI-LNSAWOW TWOTLAYOSHL €-v JYNIIS 


( 2-OTX) HONI/SNYETObY-SuNLBAYNI 


00°O 


Sroma| re 


SdIM*L4I-LNAWOW 


' 
2! igor 


——_ = > a) es 


--— 


i mentee a a ee (-—— amet eee al meee 


st a) 
\ 


a es 
matt oe 


CS-htki HO SL AC. Sey ye ry 
LY api A? Pes! NOP. oor: asst eT ale 


_ - - _ 


| “gayeae ScRSi 20% seaentd sautageen- Tame, JAuLTaRGRAT . 8-8 anu - 
_ (MRSAAG MiahTe-eeBee GST AAO!) = 
a ee a _ 7 ee a 
_ _ 


43. 


S08) USeac 


Seon ce! 


(WYYDVIG NIVULS-SS3ULS 03Z11v40I) 
SATUAS S°9TE2T YOS WVYDVIC SUNLVAUNI-LNAWOW WIILIYOSHL $-b JUNI 


( 2-01X) HONI/SNBTOYH-AUNLYAUND 


se 0 


0c*0O 


8) 


01°0O 


“S070 


ANS 


98°T=°v 


00 °S¢ 


00°0S 
SdIM° LI-LNSWOW 


OCR Sr 


00°00! 


OO Se¢l 


_ 
ie 
i 


- yeh _ 


i 
‘ ‘ : ~~ 
a —— nn Gee a ; 
—— a and aa ——— —_— —s a=58 : a > 
= = ys 7 . > 2 | itz ahs ge? oe -- oles : 
€ - ira, iT =a 7 - 


a be 2.7% 0") | 1, ons ee ; | I Tus - 7 
\ Ad '4apt } a4 se ® + ti tb a i 4 


: : - a a We. Le =) 0 oP ioe - 
y mF r aeetit ae Pe Gat (eo DS abe IGS es SESS TARGSNT dab @ wida 
we wPadGsad - ot ee. \ap7 ear? 7 . 
(MASUAIC ot Weicocas aT a }: = 


te ; a a ne ae 


(WWYOVIG NIVULS-SS3ULS G3ISIGOW) 
SAIYSS 9EAMCZT YOS WWYOVIC SYNLVAYNI-LNSWOW TWOTLSYOSHL =S-p Jundta 


44, 


( 2-O1X) HONI/SNUTO8Y-SAYNLEAYND 
SE‘O oso 80 02°0 St°0 O10 so‘o(00°0, 


a ee 
00 


elea) ie)! 


00°02t 00°09 
SdIM" L4-LNAWOW 


00°O8t 


00 ‘one 


00°00€ 


+9 aa 
2.) Ga awe 


? 
oy 
Ae 
> 
es « 
‘& 
= 
=. 
43 
oy 


45. 


(WWHOWIG NIVULS-SSIYLS GAIsICOW) 
SIINIS G°OLAZI YO WVYOVIG SYNLVAYNO-LNSWOW TWIITLAYOSHL 9-b AYNOIS 


( ¢-O1X) HONI/SNYTObY-SYNLYBAYND 


02°0 09°0 0s‘O Oho 0g "0 i oe.0 OL°o 00°O_, 
: : 
/ (@) 
S°9o gel 
ny 
pe 
(es! 
air < 
=) 
ie 
m 
PL, 
Shi 
Oat 
i 
= 
“Uy 
~~ 
S < 
0°0= S 


20st= 4 


Oe 


}——}——-}--—1 
00°001 


00°Set 


a 
= 7 


ar ® 
_ : 
a 


ie 


: 1. ee 
[ASt EL@ERS! AUT PANBAIC BRITS IAS. TGMCA SOIT IO. 
22k (WSQATO KERRI S-CORNe GRSTGGR, 


46, 


S0at 


Soahiet 


(WYYOVIG NIVYLS-SS3YLS GIISICOW) 
SATYSS 9EAMCT YOS WVYOVIC NOILVLOU-LNAWOW TIWOILSYOSHL Z2-b SYNOIS 


( ;-O1X) SNUTOUY-NOTLYLOY 
06°0 SL°0 09°0 sh°0 0€°O S1°0 


00°08t 00°021 00°09 
SdIYLS-LNSWOW 


00 *Ohe 


00°00€ 


307 SORIA OETA TOS Tem kp eaRnaiN 


6 2 @ =H af cared. 
De on tke | 4 Rage a 


47. 


(WWYOVIG NIVYLS-SS3IULS GIIIIGON) 
SSTYAS S°9TE2T YOS WVYOVIG NOILVLOU-LNAWOW TWOTLAYOSHL 8-b JYNDIS 


( t-O1X) SNBTOBY-NOT LY LOY 


IZ70 0 Bt STO 210 600 90°0 €0'0 00°C. 
TO 
Lo ] 
S°*9STG2l 
[av] 
sd 
eS 
(oun) 
on 
ro 
[om] 
S) 
~J 
i) 
0°0=5Y 9 
AP ia = 
= 
98° T=°¥ 8 


30°Set 


SdIM* L4I-LNAWOW 


48, 


(WWHOVIG NIVYLS-SSIULS G3I4IGON) 
SATYAS 9EAMCT YOS WYYOVIG NOILI91430-dVO1 WOILSYOSHL 6-% 3YNDIS 


SSHINI-NOTLIS1450 SNI1 3YLN39 
OS"€ 00°€ 0S*2 00°2 OS*} 00°1 0S°0 00°0 


cee 


00°06 00°09 00 °0€ 
SdIM-Q801 GSI Idd 


00°Oel 


00°OST 


' eS 


. 2 
Et) 
> 
- 
c= 
™ 
ww 
>? 
a) 
yh 
~~ 
ye 
> 
we 
OS 
hy 
. 
aay 
oa 


es 3 


SPORE ge EDSELIceT | 


ne — 


49, 


(WWYOVIG NIVYLS-SS3YLS GIISICOW) 
SdIUaS G°9TE2T YOS WVYOVIC NOILOATSS0-dVO1 WOILSYOSHL Ol-b JYNDIS 


SSHINI-NOILIS1S30 SNI7 3YLN39 


02°0 09°0 0S°0 Oh°O 0€°O 02°O OL°o 00°O_ 
+= 
SS 
S°9Tgel 
on 
a5 
-=0 
ae 
— 
m 
ea 
(Sb) 
Sir 
QaQ 
OD 
i. 
=< 
S a3) 
0°0= V Hu) 
(en) 
[en] 
S 
COST =. 
S (op) 
el = y = 
fom] 
(an) 


00°SL 


: “» + @ PAs 
Ti7,{ CA 


4 
c@.) 
wel 
=~. 


a) > 
ae rar) =: 


— a 
i. 4 , 
=~ 
% 
ae > 
2 
Qual 
haa 
«4 
2 
= 
t 
ye 
7 * 
par i vs 2 
fi of a 
= / ww rl 
- , 
: " ZA , a . 
- P 


Pi ft 
d , ol 
a 
a. 
_ 


_ 


_ . Bier 


yes 
) 
‘Lt 
a0 


Nelo | a 
7 rm bt . 


ad 


90. 


(WYYDVIG NIVYLS-SS3IYLS CG3I4IGOW) 
SATYuaS 9EAM2T YOS WVYOVIC NOILVYIOT SIXV IWYLAAN-LNSWOW TWOILSYOSHL = TI-p ayNdIs 


SdIM*LI-LNSWOW 
00 *00€ 00 "one 00°081 00021 00°09 00°0 


00°O 


cedMed 


00*l 


SSHONI-NOILYDO)] SIX IWYLASN 


8T°T = V 


00*e 


00° CcO°e 


00°S 


teres 
tte i. Sell VAP 
a : 

s s@ %. 


Le 


43 


iHEOHEL cw ven ll tie 


nee aay 
= : 
ae 


- 
¢ 
cana? 
aii ww = 
———— 


Tice 


~ 


4) 


, 


aha ‘ © > 
tit} JEWL- et 


_ 
ae) 


ye : L2Pc s” § ; 
re pn et cr Na i 


= 


ai. 


51. 


(WWYOVIG NIVULS-SS3NLS GII4IGONW) 
SSTYSS G°9TE2T YOs WWYSVIC NOTLVYIOT SIXVY TWWYLNSN-LNSWOW TWOILSYOSHL 2I-b 3yNdI4 


SdIM*LI-LN3AWOW 
00°S21 00°00! 00°SZ 00°0S 00°S2 00°0 


00°O 


Seale 


00 °e 00° 


SSHONI-NOTLY9O1 SIX¥ IWYLNAN 


00°€ 


CO*h 


00°S 


cae~ previa Cr verue! 


(WOOTRIED are 


a a 
toy ers 


a) wiht 


nO 


od 
oad 
a 
. 7 
Fe 
> 
&« 
& 
hae 
wat 
a 
‘s - 
hy -* 
= 


aie eee “ 


Taal 


A 


52. 


I2B1 6552 SERTES 


12WF36 SERIES 


T2ViPCT SERIES 
12WF31 SERIES 


0 50 100 150 200 


A. Fis (Kips) 


FIGURE 4-13 EXPERIMENTAL RELATIONSHIP BETWEEN MOMENT RATIO Mare ait AND 
YIELD FORCE IN SLAB REINFORCEMENT A, Fis 
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FIGURE 4-14 THEORETICAL RELATIONSHIP BETWEEN MOMENT RATIO Myy¢/My7¢ AND 
YIELD FORCE IN SLAB REINFORCEMENT A, Fy. 
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CHAPTER V 
DISCUSSION OF RESULTS 


SA INTRODUCTION 

The following discussion centres on the results of the 
analyses based on idealized and modified stress-strain disgrams. 
Comparisons are made between theoretical and experimental results. 
Factors which may have affected both theoretical and experimental 


results are discussed. 


Oe COMPARISON OF TWO BASES FOR THEORETICAL SOLUTIONS 

The flexural behaviour of a composite member depends on 
cross-section properties, material properties and the type of loading 
applied. Material properties may be obtained from a tension test on 
a coupon of the material. The results of such a test may be defined 
in the form of an idealized stress-strain diagram. The idealized 
stress-strain relationship, however, does not satisfactorily define 
the response of the member to applied load. . The presence of residual 
stresses, in the steel section as a result of the rolling process, 
is a factor in the resulting stress distribution due to applied 
moments. The residual stress pattern in the flange varies from 
compression at the flange tips to tension at the flange-web junction 
A coupon test does not provide information about residual stresses, 


since upon cutting the coupon these stresses are released. The 
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57} 
modified stress-strain diagram attempts to approximate the effects 
of residual stresses. 

'For the theoretical analysis, material properties such as 
yield stress, yield strain, strain at strain-hardening and strain- 
hardening modulus were obtained from coupon tests which were made 
in connection with the experimental investigations of composite beams 
at the University of Alberta. 

The general flexural behaviour: of a composite beam may be 
defined by moment-curvature relationship. The theoretical moment- 
curvature relationships (FIGURES 4-3 to 4-6) are very similar in 
nature to those obtained experimentally (FIGURE 4-1 and 4-2). Use 
of the modified stress-strain diagram in the analysis results in 
curves which bear a closer resemblance to the experimental curves 
than those obtained with the idealized stress-strain diagram. 
Generally, the modified stress-strain diagram and the idealized 
stress-strain diagram result in almost identical vance of 
ultimate moment. The curvatures at ultimate, however, are significantly 
different as the plastic plateau is neglected in the modified stress- 
strain diagram. The unloading portion of the moment-curvature 
relationship, subsequent to ultimate moment, is not defined by the 
theoretical analysis. The effect of increasing the amount of 
longitudinal slab reinforcement is similar for both experimental 
and theoretical relationships. Relationships derived from theoretical 
moment-curvature relations appear in FIGURES 4-7 to 4-12. 

In order to fairly assess the validity of the theoretical 


analysis, one should first examine both the theoretical and 
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58, 
experimental values and decide if any of them are of a questionable 
nature. TABLE 4.1 indicates that BEAM 25 does not follow the trend 
of other beams in the same series with respect to ultimate moment 
value. Therefore, the experimental results for BEAM 25 are suspect. 
Similarly, the 12WF27 series should probably be excluded from the 
comparison since the material properties defined by coupon tests 
appear to be in error, In the 12WF36 series BEAM 18 is inconsistent 
with respect to the deflection at ultimate moment when compared 
with other beams of the same series. During testing, BEAM 18 
developed a twisting rotation and, therefore, had to be unloaded so 
shims could be placed under the loading bridge before the load was 
reapplied. Ultimate deflections rather than ultimate curvatures 
are compared as it is felt that experimental values of deflections 
are likely to be more accurate than values of curvatures. At moment 
values above Mj, strain gauges became inoperative. Therefore, 
experimental curvatures at moment values above M, are based on strain 
gauges located close to the neutral axis. 

The theoretical values of ultimate moment differ from 
experimental values by a maximum of -8% for beams in which the 
compression flange strain-hardens prior to failure by local buckling. 
For sections which do not reach strain-hardening of the compression 
flange, that is, the composite 12B16.5 series, the theoretical 
ultimate moment value is 12% below the experimental value. For the 
12WF31 series, deflections based on the idealized stress-strain 


diagram are up to 80% greater than the experimental deflections 
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whereas the values based on the modified stress-strain diagram show 
a difference from experimental values of -16% to +14%. The 12WF36 
Series show a maximum difference of -26% between experimental and = 
theoretical ultimate deflections based on the idealized stress- 
strain diagram. Based on the modified stress-strain diagram the 
maximum variation is -50%. For the 12B16.5 series, use of the 
idealized stress-strain diagram results in a 47% difference while 
the modified stress-strain diagram gives a deflection 53% under 

the experimental value. The above comparisons indicate that the 
value of ES or Bea affects the ultimate values of curvature, 
rotation and deflection to a greater extent than it affects the 
ultimate moment value. This is also observed in results obtained 


with ES or Ee Vanvind TrOMmmaEALOO tose /cOroratronmse< 102051" 


t 
to 10 x 10°psi. From FIGURES 4-3 to 4-6 it is also evident that 

a small increase in the moment or load value, for values above M, 
results in a significant increase in curvature and related values. 
Based on theoretical results for the 12WF31 series it appears that 
the modified stress-strain diagram results in a better agreement 
between theoretical and experimental values at ultimate moment 
than does the idealized stress-strain diagram. 

On the basis of all the comparisons, the modified stress- 


strain diagram, therefore, appears to provide the better solution 


for the flexural behaviour of composite beams in negative bending. 


Died MATERIAL PROPERTIES 


The stress-strain diagrams are derived from values of yield 
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60. 
stress and strain, strain-hardening strains and strain-hardening : 
modulus obtained from tension tests on coupons taken from the beams 
tested in the experimental investigations. As TABLE 4.6 shows, a 
wide variation can exist in the above values obtained from coupons 
taken from the section. It is, therefore, questionable whether 
average values obtained from coupon tests are likely to be more 
valid than standard values specified for the given material, in this 
case.G40.12 steel. The strain-hardening modulus, in particular, 
appears to vary considerably. Other investigators indicate that 
the strain-hardening modulus may vary between E/100 and E/30, in other 
words fromse<110rosi. to 10 x 10°psi. 

Theoretical results based on standard material properties 


for CSA G40.12 steel (F, = 44 ksi, E = 29000 ksi) with Bee and Ect 


y 
ranging from E/100 to E/30 were obtained. The results indicate that, 

as Et is increased, ultimate deflections decrease whereas the moment 

at ultimate remains fairly constant. Compared with experimental results, 
the moment at ultimate is over estimated by a maximum of 15% while 

the deflections at ultimate represent upper and lower limits to the 
experimental values. Based on the above, it appears that standard 


values may be used to determine the flexural behaviour of composite 


beams in negative bending. 


5.4 NEUTRAL AXIS LOCATION 
FIGURES 4-11 and 4-12 indicate the initial tendency of the 
neutral axis to move towards the tension flange with increasing moment 


and increasing longitudinal slab reinforcement. Initially the neutral 
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Gl 
axis is located at the centroid of the composite section. Its 
location is constant until the moment reaches the value Mo. -Once 
yielding begins, the neutral axis moves towards the tension flange 
with increasing moment. This movement is halted either by strain- 
hardening of the compression flange or yielding of the longitudinal 
Slab reinforcement. The latter is more likely to happen before the 
former. As moment increases still further, the neutral axis moves 
towards the compression flange. For the 12B16.5 series, which do not 
exhibit evidence of strain-hardening, there is no tendency for the 
neutral axis location to shift towards the compression flange at 


moment values above Mo: 


Berd ULTIMATE LOAD CONDITIONS 

Values in TABLES 4.1, 4.2 and 4.3 have been used in several 
plots to show the effects of longitudinal slab reinforcement and size 
of the steel section on the ultimate moment capacity. In FIGURE 4-14, 
for the 12B16.5 series it is evident that an increase in the amount 
of longitudinal slab reinforcement results in a substantial increase 
in the Sar mae 
gradual increase in Miae/Mare for low amounts of longitudinal slab 


ratio. The other series, however, exhibit a more 


reinforcement. Above the limit of reinforcement which places the 
neutral axis into the tension flange, the increase in Muie/Mort is 

1% or less as the area of reinforcement is reduced to the upper 

limit value given by EQUATION (3-1). The increase in the experimental] 
Muie/MGi¢ ratio is less than Lie ae increase in the Mij4/Mgi¢ ratio 
is not as great as the increase in the Mo/ Mo ratio for the series 


involving wide flange sections. The increase in ultimate moment is, 
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62s 
therefore, not directly proportional to the increase in the simple 
plastic moment (as shown in FIGURE 4-15). FIGURE 4-13 shows the 
experimental results for comparison. 

FIGURE 4.16 indicates that the theoretical moment capacity 
is dependent on the b/t ratio of the steel section. For composite 
sections with similar A./ Awe ratios, where the compression flange 
strain-hardened prior to ultimate moment conditions, an almost 
linear relationship exists between the b/t ratio and the ultimate 
moment. This relationship is verified by the experimental results. 
BEAM 42 does not follow this trend since local buckling occurs 


prior to strain-hardening of the compression flange. 
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CHAPTER VI 
CONCLUSIONS 


From a comparison of theoretical and experimental results 

it is concluded that: 

1) The ultimate moment capacity of a composite section in 
negative bending may be computed by the proposed theory. 
Theoretical values of ultimate moment, based on actual material 
properties as obtained for test specimens at the University 
of Alberta, are a maximum of 8% below experimental values 
for the beams failing after strain-hardening of the 
compression flange. For section failing prior to strain- 
hardening of the compression flange, the theoretical 
ultimate moment is 12% below experimental values. 

2) Use of a modified stress-strain diagram consisting of an 
elastic range and a Sea ae eemnne range produces results 
closer to test results than those obtained using an idealized 
stress-strain diagram consisting of elastic, plastic and 
strain-hardening ranges. 

3) The upper limit on the effective area of longitudinal 
slab reinforcement can be satisfactorily determined on the 
basis of a neutral axis location, for the plastic moment 
value, at the inner face of the tension flange. 


4) Ultimate moment values are not particularly sensitive 
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to the value of the strain-hardening modulus, but curvatures, 
deflections and rotations are very sensitive. Theoretical 
results obtained using the modified stress-strain diagram 
with the strain-hardening modulus varying from E/100 to E/30 
result in a relatively small increase (less than 1%) in the 
ultimate moment whereas corresponding curvatures decreased 
from an upper limit for the experimental values to a lower 
Tame. 
The theoretical analysis verifies the following conclusions 
established in experimental investigations by Lever (3) and 
Davison (2): 
1) For a given steel section, significant increases in the 
ultimate capacity of composite beams in negative bending can 
be achieved by the addition of longitudinal slab reinforce- 
ment. 
2) The increase in the ultimate moment capacity is not 
directly proportional to the increase in the simple plastic 
moment value. 
3) For a given steel section, an increase in the amount of 
longitudinal slab reinforcement results in a reduction in 


the rotation capacity of a negative plastic hinge. 
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APPENDIX A 
COMPUTER PROGRAM 
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A.l COMPUTER PROGRAM 

The computer program, written in FORTRAN IV language for 
the IBM 360/67 system, calculates and plots moment-curvature, 
moment-total rotation and load-deflection relationships and 
determines the location of the neutral axis for composite sections 
in negative bending. The calculations follow the procedure described 
INPCHARTERSIII.< 

Subroutines are used to calculate the centroid of the 
section and the moment of inertia about the x-axis, to list most of 
the calculated values and to plot the calculated relationships. The 
flow chart for the main program is given in APPENDIX A.3 and some of 
the notation used is described in APPENDIX A.2. The formulation of 
the subroutines can be followed in the program listing presented in 


APPENDIX A.4. 
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Hee NOTATION USED IN COMPUTER PROGRAM 
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SPRY Reo TRY We oTRYVR 


SolbseolsiW,, olotk 


LPEOT 


FACTOR 


DIST 


XDIST 


GE 


steel section and physical 
properties of steel section 

brief description of input 

yield stresses of web, flange and 
longitudinal slab reinforcement 
Young's modulus for flange, web 

and longitudinal slab reinforcement 
strain hardening modulus for flange, 
web and longitudinal slab reinforce- 
ment 

yield strain of flange, web and slab 
reinforcement 

strain hardening strains of flange, 
web and slab reinforcement 
determines scale to which plots are 
drawn in Subroutine GRAPH 

amount of slip considered to occur 
between slab and steel section 
distance between reaction and load 
point 

length of beam, this is divided by 

2 and value obtained is also 
designated as XDIST 

elastic shear modulus 


Poisson's ratio 
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A-4 
AF, AW - ratio of Young's modulus to strain 


hardening modulus for flange and web 


GTE, GIW- - inelastic shear moduli of flange 
and web 
WI - variable used in calculating length 


of flange required to yield for 


formation of local buckle 


ICHECK - integer counter 

AREA 1, AREA 2, AREA 3, - variables used in calculating 

AREA 4, AREA 5 rotation and deflection of beam 

Fog Fry FW - forces in reinforcement, flange and 


web used in computing simple plastic 
moment 
NA - neutral axis location, above 


compression flange, for simple plastic 


moment 
¥ - depth of top flange in compression 
MSS MTF, MBF, MW - contributions to plastic moment of 


reinforcement, top flange, bottom 
flange and web 

MP - simple plastic moment 

Zz - assumed distance between neutral axis 
and compression flange 

INCR - value by which neutral axis is 
shifted if assumed location of 


neutral axis does not result in 
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FCW, FCYW, FCEW 
FCF 
FTW, FTYW, FTEW 
FTF 
FAS 


MUAS, MUTF, MUTW, MUTYW, 
MUTEW, MUCEW, MUCYW, MUCW, 
MUCF 

MULT 

AMPS 


APHI 


XOPT 


XMPS 


A-5 
equilibrium of section 
compressive strains in web and flange 
tensile strains in web, flange and 
reinforcement 
elastic portion of web in compression 
depth of compression zone of web which 
1S inelastic 
forces in web in compression 
force in compression flange 
forces in tensile portion of web 
force in tension flange 
tensile force in reinforcement 
sum of forces in section 
integer counter 
contributions to moment in section 


due to specified curvature 


total resisting moment of section 
moment at first yield of compression 
flange 

curvature at first yield of compression 
flange 

distance required for local buckle 

to form 

ultimate moment (moment at which 


local buckling of compression flange 
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SK 


ADIST 


PINCR 


JJ, Kd 
SRATIO 
BMTOT (JJ), BPHI (Kd), 
ROT (JJ), ALOAD (JJ), 


DELTA (JJ), ADIST (JJ) 


PHI, UNIT 


XRATIO 
AIX 


YBAR 
NNKK 


occurs ) 

torsional spring constant 
distance between neutral axis and 
centroid of steel section 

amount by which curvature is 
increased 

integer counters 

inelastic portion of beam 

related values, at any loading, 
of moment, curvature, rotation, 
load, deflection and distance 
between neutral axis and centroid 


of steel section 
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curvature of section at a specified 


load 

elastic length of beam 

moment of inertia of section 
considered 

centroid of section considered 


integer counter 
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FLOWCHART FOR MAIN PROGRAM 


READ SECTION AND 
SECTION PROPERTIES 


WRITE SECTION 


HE 


NO 
READS TE 


READ YIELD STRESSES 
OF WEB, FLANGE AND 
SLAB REINFORCEMENT 


niceties Das oc aa) a5 


READ ELASTIC MODULI 
AND STRAIN HARDENING 
MODULI 


CALCULATE YIELD 
STRAINS OF WEB AND 
FLANGE 


READ STRAIN HARDENING 
STRAINS OF WEB AND 
FLANGE 


READ MATERIAL 
PROPERTIES OF 
LONGITUDINAL SLAB 
REINFORCEMENT 
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READ LPLOT, FACTOR, 
BBS ere GBR 


NOLS T= DIS i 


GE = 11200 
CALCULATE AF, AW 
US=sUea 


CALCULATE GTF, GTW 


READ AREA OF 
LONGITUDINAL SLAB 
REINFORCEMENT 


PREPARE SOME OF THE 
OUTPUT, ICHECK = 1 
CALL CALCI 

CALCULATE SIMPLE 
PLASTIC MOMENT 

WRITE SIMPLE PLASTIC 
MOMENT 


CALCULATE EFFECTIVE 
YES-MAREA OF SLAB 
REINFORCEMENT 
NO 
INITIALIZE SOME H 
VARIABLES 


CALCULATE YIELDED 
LENGTH REQUIRED FOR 
LOCAL BUCKLE 
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INITIALIZE UNIT 


DETERMINE STRAIN 


DISTRIBUTION OVER 
CROSS-SECTION 


COMPUTE STRESS 
DISTRIBUTION AND 
CORRESPONDING FORCES 
DETERMINE SUM OF 
FORCES 


NO 


INCREMENT NEUTRAL 
AXIS 


Ke 


CALCULATE MOMENT 
OF SECTION 


CALCULATE XMPS 


7 
Th 
T 


44 


WRITE XMPS, XOPT, Z, 
MULT 


CHECK eae 


INITIALIZE SOME OF 
THE VARIABLES 


PHl) = PHI ?seRINCR 
UNIT = PHI 


|smToT (Jd) - XMPS| vesi( 700) 


<0.05 


(Chg 1 ES) 
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A.4 PROGRAM LISTING 


OIMENRSTICN DATA(C1024) 

CIMENSTON SRATIO(4000) 

CCM¥MCN KOT (4000) 

CCOMMOGR JPLOT »LPAGF »eTITLE( 20) eLPLCTsSECT(10)sXMPS,ADIST( 4000) 
CCMMON AC10)+T(10) owW010)5E019)2ESH(10) »S(10)+.SSH(10) -SNTOT( 4000) 
CCMNCK BPHI(4C00) »sJJeKJeAL0AD( 49000) eDELTA( 4000) eAIX 

REAL MUTF s4UTY &sMUAS oMUTW eMUCF eMUCY We MUCWeMULT 2s INCRsMUCEWe MUTEW 
REAL MSs MTFeMBFoMWoMP MTF TseMTFCoMWT eMWCoNA 

CALL FLOTS (DATA(1)24096) 


JPLOT=1 

Cc *% 
Cc ** READ IN SFCTICN AND SECTICN FRCPERTIES 
Cc ** 
Cc xe 
Cc ** D=DEPTH OF X-SECTICN, B=wIOTK OF FLANGE» tTT= THICKNESS GF FLANGE® 
Cc ** Whe= THICKKESS OF WEB 
E ** CO FQUALS CISTANCE THAT REINFe OF SLAB IS ABOVE TENSION FLANGE 
(a xe 
GC ae 

1 PFEAC(S2100) (SECTCI)»T=122)e0DeBeTT»wWeDD 


100 FCRMAT(2X%s2A4%55F10 06) 
WRITE( 66200) (SECT(I) eT =122) 
200 FCRMAT(1IH1Le¢10Xe*SECTICN IS %22A4e//) 
IF(DeLE°0-0)GU TO 909 
. READE, 623)(TITLEC IL) »l=1+20) 
623 FCRMAT(20A4) 
WRITE( 6,624) (TITLECI) sl =1220) 
624 FORMAT(C12X%s20A4e/4/) 
eg 
** READ IN YIELO STRESSES OF WEE.sFLANGE AND REINFORCEMENT 
** TY = YIELD STRESS OF WEB» TYF = YIELD STRESS OF FLANGE 
** TYS = YIELD STRESS GCF REINFORCEMENT STEEL *¥**¥4% 
** 
REAL (5.300) TYWseTtYFoTYS 
300 FCRMAT(3F10-6) 
H=0-(2e*TT) 
( ** 
E ** READ IN ELASTIC MODULUS AND STRAIN HARDENING MODULUS 
(e oe 
REAC(5:500) EF sEWeESTFeESTW 
500 FCFMAT(4F20-6) 
G ae 
G ¥%# CALCULATE YIELD STRAINS OF WEB AND FLANGE. 
G 4 
STRYF=TYF/EF 
STRYW=TYW/EW 
Cc *% 
(e **# READ IN STRAIN KARDENING STRAINS 
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A=13 


G ** 
REAC(£+4000)STSIF »STSTW 
4000 FCRMAT (2F20.10) 


READ IN MATERIAL PRCPERTIES FOR SLAB STEEL 


READ(ES +4002) ERsESTReSTSTR 
4002 FORMAT(2F20.210) 
STRYR=TYS/ZER 


LPLOT OETERMINES TRE SCALE TO WHICH THE MOMENT ROTATION CURVES ARE DRAWN 
FACTCR EQUALS SLIP BETWEEN REINFeAND STEEL SECTION 
OIST EQUALS DISTANCE BETWe REACTICN POINT ANO LCAD POINT 


(ay tal tay tay (a) 


REAC(££.4001)LFLCT»FACTCR»DISTeXOIST 
4001 FCRMAT(T2s3F10c¢&) 
XCIST=XOIST/2a4 


GCE=112006 
Cc ** . 
Cc 4% 
(E “ee AF = RATIO CF ELASTIC MODULUS CF FLG TO STRAIN HARDENING MODULUS 
Cc ** AW = RATIO OF ELASTIC MODULUS OF WEB TO STRAIN HARDENING MUDULUS 
Ce ** 
(e ** 
AF=FF/ESTF 
Aw=EW/ESTW 
(e *** 
G *% U = POTSSCNE RATIO 
G ** 
U=e3 
fe ** 
(a ae 
(e ** GTW, GIF ARE INELASTIC SKEAR MCDULTI 
GC *% 
Cc xs 
GTIW=(2eXGE)S (Let Aw/(4e%*(1.e+U)) ) 
GTFH(2Ze*GE)S (1 e tAFS(4e¥*(12-+U))) 
WI=CC(E#*3)XTT #42) F70)S( 160 #1440) 
G ** 
Cc ** READ IN AREA CF REINFORCEMENT 
G ** 


11 READ(E e400) AS 

400 FCRMAT(F1I0.6) 
AREA1=O6 
AREA2Z=06 
AREA3=06 
AREA4=006 
AREAS=06 
ICKECK=1 
ITF CAS2GTe1020) CO TO 1 
DO 80C K=1+5 
IF (KeEGe1)GO TO 201 
TECK eEGe2)GO TC 802 
IF (kK eEQe3)GO TO 803 
IF (KeECe4)GO TO 802 
TECK el Ges) GU TTC e048 

801 CCNTINUE 
A(K)=4S 
T(K)=00*2e 
w(K )=2e 


WA 
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802 


£03 


804 


E(K)=ER 
ESH(K)=ESTR 
S(K)=STRYR 
SSH(K)=STSTR 
GCC TO 805 
CONTINUE 
ACK) =E*TT 
WK =—T Fi 
w(K)=8 
E(K)=EF 
ESH(K)=ESTE 
S(K )=STRYF 

© SH(K)=STSTF 
GC TQ 805 
CCNTINUE 
A(K)=H* ww 
T(K =F 
W(K)=)W 
E(K)=FW 
ESt(K )=ESTW 
S(K)=STRYW 
SSH(K)=STSTW 
GES TOsSsoS 
CCNTINUVE 
A(K)=0¢ 
T(K)=064 
W(K)=06 
E(K)=06 
ESH(K)=0Oe 
S(K)=Ce 
SSH(K )=06¢ 
CCATINUE 
CONTINUE 
WERITE(E.919) 


FCRMAT(10X, '-~-~~-------~------- ----------------- --------------- *) 


CALL CALCTI 

WRITE (€@s600)AS 
FORMAT(10Xe?# AREA OF REINFORCEMENT IS 
a% 

** COMPUTE PLASTIC MCMENT ¥**% 
xe 

FS=AS4TYS 

FF=TYF*B¥TT 
FWw=TY**HeWw 
IF(FSeGTeFW) GO TO 2 
IF(FSeLTe-FW) GC TO 3 
MS=FS*(TT+00) 
MTF=FFR(TT/2 0) 
MBF=FFeCHETT/20) 
Ma=FWe(H/20) 
ME=(MS4tMTFtMBFtMW)S 126 
NA=H 

GC TO 4 

CCNTINUE 
F=(Fu4+(2e*FF)-FS)/26 
Y=(F)/(TYF*B) 
FFTC=FF-F 

MS=FS*(O0DtY) 
MTIFT=F*(Y/20) 
MTFC=FFTC#((TT-Y)/20) 
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SQUARE 
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MMW=FW#(H/204+(TT-Y)) 
MBF=FFX(H+( (CTT /2¢)+4+(TT-Y))D) 
MP=(MS+MTET+MTFC+MWtMBE D/120 

NA=Ht(TT-Y) 

GC TO 4 i 
CCNTINUE 
F=(FW-FS)/2e 

Y=(F)/S(TYW*WW) 

FYXC=FR—-F 

MS=FS*(Y+TT+0D) 

MYIF=FF*CY+(TT/20)) 

MEF=FFX*(H-Y#+(TT/72e)) 

MW¥T=F*(Y/20) 

MwWwC=FRCK((CH-Y)/20) 

NF=(MS4+MTF4MBE tMUT4+MWC S120 

NA=H-Y 

CCNTINUE 

WRITE (65501 )MP 


FCRMAT(L2X,* PLASTIC MOMENT OF SECTIUCN IS "5 F8e30"FT eKe®) 


WRITE (Es700)NA 
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FCRMAT(12Xs* NEUTRAL AXIS LOCATICN IS #8 »F7e3s* INCHES ABOVE COMPRE 


SSCiCN FLANGE ® s/) 
IF (NA eGTeH)AS=FWU/STYS 
WRITE (69297) AS 


FORMAT(12X-e*AREA CF STEEL RECUIRED [IS ®,F6235* SQeINe®) 


Z=NA-1e 

TF (CAS eEQ2020)Z=t+l2e 
TNCR=e2 

ELIST=0O6 

SULM=0e 

N=0 

CCNTINUE 
SK=(GTFeWW**3)/7(03 242) 
XL=3-214159*SQRT(SQRT(ESTF#WI/SK) ) 
XxCFT=2 e*XL 

UNI T=STRYF/(ZtTT/20) 
CCANTINUE 
ELAST=STRYW/SUNIT 


DETERMINE STRAIN OISTRI@GUTION OVER X-SECTION 


ECF=UNITH(Z+TT/20) 
ECW=UNIT*Z 

ETW=UNIT*(H-Z) 
ETF=UNITX(H-Z4(TT720)) 
EAS=UNIT*(H-Z+TT+C0) 4FACTOR 
IF(ECW.LE*STSTW)GO TO 60 


ZZ = CEPTH OF CCMPRESSICN ZONE OF WEB WHICH IS PLASTIC 
Z = CEPTH OF WEE IN CCMPRESSICN 

BUASiIG = EEAST IG REGTCN ER Wes 

CCNSICER COMPRESSICN ZCNE OF WEB 


Z2Z2=((STSTW-STRYW)¥(Z-ELAST))/7(EC#-STRYW) 
TCW=(ECW-STSTW)¥ESTW 

FCw=(TC www (Z-ZZ-ELAST) 7260 

ele, fer fairl 

CCKTIRUE 

IF(EFCWeLE-STRYW)GO TO 601 
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601 


602 


62 
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€3 


€4 


(eS) 


641 


€42 


€6 
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FCW=06 

CCATIAUE 
FCYW=TYW*WW*(Z-ELAST) 
FCEW=(CTYW*¥WW*ELAST)/26 
GC TO 602 

CCNTINUE 

ELAST=Z¥1le 

FCw=06 

FCYW=06 
FCEW=(ECWHEW*Z WW) 20 
CCNT INUE 


CCNSIOER COMPRESSICN FLANGE 


IF (ECF eLEeSTRYFICO TO 92 
IF (ECFe«LE*«STSTF)GO TO 62 
TCF=(ECF-STSTF)*XESTFATYF 
FCF=TCF*B*TT 

GC TO €3 

CCNTINUE 

FCF=TYF*B*TT 

GC TO €3 

CCATINUE 

FCF=B*TT#ECFE*EF 

CONTINUE 


CCNSICER PORTION OF wEB IN TENSICN 


IF(CETwWeLEeSTSTW)EGO TO 64 
TTw=(ETw-STSTh)*ESTW 
FIw=(TTW*Www* (b}-Z-ZZ-ELAST) ) 720 
GEC TO 65 

CCNTINUE 

IF(CET¥ eAEoSTRYW)GC TO E41 
FT™Ww=0Oe 

CCNTINUE 
FIYs=TYw*Wwt(H-Z-ELAST ) 
FIEW=(TYWtWw*ELAST)/26 

GC T0 &€42 

CCNTINUE 

ELIST=H-Z 

FTw=06¢ 

FTYW=06 
FTIEW=(ETW*eWWeEWw*ELIST)/20 
CCNTINUE 


CCNSICER FLANGE IN TENSICN 


IF (ETF LEeSTRYF)IGO TO S3 
IF CETFeLE*STSTFIGO TO 66 
TTF=HC((CETF-STSTF) *¥ESTF ) tT VF 
FIF=TTF*B*TT 

GC TO 67 

CCNTINUE 

FIF=TYF*B*TT 

GC TO 67 

CCNTINUE 

FTF=O*TT *eETF EF 

CCNTINUE 
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CCNSTIDER SLAB REINFORCEMENT 


IF (FAS eLE2-STRYF)CGO TO 68 
FAS=TYS*AS 

GC TO 69 

CCNTINUE 

TAS=EAS*ER 

FAS=TAS*AS 

CCNTINUE 


TEST FOR EQUILIBRIUM OF X-SECTICN FORCES 


SUMF=(FAStFTFLFTwtFTYwt+F TEW )-FCEW-FCY w-FCW-FCF 
N=N+t+1 

IF CACS(CSUMF) eLTee1)GO TO 17 

IF(NeLT22)GO TO 79 

IFCSUMF.LT20-0)GO TC 70 

IF ( (SUM—-SUMF ) eLT e0e0eANDoSUMeGT.2020)60 TO 71 
IF ( (SUM-SUME ) «GT e0e0eANDeSUM.LGT2020)GO TO 73 
GGT ems 

CCNTIAUE 

IF ( (SUN—SUMF ) eLT e0 eOeAND2SUMeLT2020)GO TO 75 
IF (( SUM—SUMF ) «GT eDeCeoANDeSUMeLTe020)GO TO 76 
GOR NO maT 

CCATIAUE /’ 

ele, Tiel 747¢ 

CCATINUE 

IFCINCReLTeO-)GO TQ 50 

INCR=-—e7*ABSCINCR) 

GC TO 19 

CCATIARUE 

INCR=e7*ABS( INCRD 

GC TO 19 

CCNTINUE 

IFC INCReLT2O006)GC TO SI 

INCR=1-S*ABSC(INCR) 

GC To 19 

CCNTINUE 

INCR=-12¢5*ACS( INCR) 

Gc TC 19 

CCNTINUE 

IFCIKCReLT-O2VGO TO S2 

INCR=1.25*ABS( INCR) 

Gc TO 19 

CONTIAUE 

INCR=-1-225*ASBS(INCR) 

Gc TO 19 

CONTIAUE 

IFCINCReLTeO-)GCO TO S3 

INCR=-1-e35*ABSCINCR)D 

Gc TO 19 

CCNTINUE 

INCR=1-35*AB8SC( INCR) 

GC TO 19 

CCNTIAUE 

IFC INCReLTeO2)GO TO SS 

INCK=-eS*ABS(INCR) 

cc Tc 19 

CCATIALE 

INCR=eS*ABS( INCR) 
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Gc TO 19 

79 CCNTIRUE 
IF(SUMF.LT-0-)GCO TC $4 
TINCR=ABSCINCR) 
cc 10 19 

54 CCNTINUE 
INCR=-12*ABS(INCR) 

19 CCNTINUE 
SUM=SUMF 
Z=Z+IKCR 
IF CICFECKeEQ-2)GO TC 777 
GC TO 763 

17 CCNTINUE 
MULAS=FAS* (H+TT+00-Z) 
MLUTF=FTF*(H-Z4+(TT/20)) 
MLTW=FTW*(( (0 (-Z-ZZ2-ELAST )*¥20¢)/30¢)+Z7Z4EL AST) 
MUTYW=F TYW* ( (H-Z-ELAST)/2-e+ELAST) 
MUTEW=(FTEW*FELAST#2 0/360 
TF(ELIST LE (H-Z))MUTEW=(FTEWXELIST#20e)/ 36 Z 
MUCEW=(FCEWXELAST*22)/30 
VUCYW=FCYW*((Z-ELAST)/26+ELAST) 
MLC w=F Cw ( (( (2-ZZ-ELAST) *20)/30)+Z2Z4+EL AST) 
MUCF=FCF*(Z+TT/20) 
MULT=(MUTF+MUTYWtMUAS+MUTW+MUCF AMUCYW tMUCWH+MUTEW+MUCEW) S126 
IFC ICRECK.-EQe2)GO TO 778 


AMFS EQUALS MCMENT AT FIRST YIELO CF COMPRESSIGN FLANGE 


ANP S=MULT 
APHI=UNIT 
XNPS=MULT*¥DIST/(OIST-XCPT) 
WRITE (6e765)XMPS 
765 FCRMAT(12Xs*ULTIMATE MCMENT BASEC ON RECALCe’ VALUE OF MPS IS *eF 8e 
VS? Fie KIePS*) 
WRITE (60773)XCPT 
773 FCFMAT(12X5*OPTIMUM LENGTH FOR LCCAL BUCKLE TO OCCUR IS *"sFT7e3e® 
1 INCHES*) 
WRITE (65774)2Z 
774 FCRMAT(12X9*NeAe LOCCATICN FOR CALCULATION OF MPS IS *sF623,° INCH 
1ES ABCVE COMPRESSION FLANGE *) 
WRITE (69775) MULT 
775 FCRMAT(12Xe*CALCULATED VALUE CF MPS IS *sF8e3e? FT.» KIPS*»s//) 
ICFECK=2 


THIS FART OF THE PRCGRAM WILL DEVELOP MOMENT-CURVATURE RELATIONS 


ADIST(1)=Z-(CH/20) 

Z=NA-1e5 

IF (CAS eL Te -GO01 )Z=HJ 26 

PINCR=-00002 

JJ=2 

K J=2 

SRATIO(1)=06 

BNTCT(1)=06« 

BPKI(1)=06 

RCT(1)=06 

ALOAD(1)=O084 

DELTA(1)=0.6 

PHI=UNIT-PINCR 
WenG CCNTINUE 
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PHI=PHI+PINCR 
UNIT=FHI 

Ge OQ Wri, 
CCNTINUE 
BMTCT(JJ)=MULT 


‘BFHI(KJ)=PHI 


ADIST EQUALS LCCATICN CF NeAe ABOVE CENTER OF STEEL SECTION 


ACIST(JJ)=Z-CH/26¢) 


DELTA EQUALS CENTERLINE DEFLFECTICK OF BEAN 


XRATIO EQUALS ELASTIC LENGTH OF EEAM 


XRATIC=CAMPS*CISTI/SBMTCT( II) 


TFCCBMTOT(CSJS)—-AMFS) eLEe e-0005)XRATIO=DIST 


SRATIC(IJI)=DIST-XRATIO 

AREA1=XRATIG* «-S*#EPEKI( JJ) 

IF (EPHI( JJ) eGT eAFHI )AREAL=KRATIC* eS *APHI 
AREA2=(XDIST-DIST)*BPHI( JJ) 
AREAZ=(SRATICO(JJ)—-SRAT ICC JIJ—1) ) * (BPHI (JS) 4+EPHI( JJ-1))*€25 
ARE AG=AREA4tAR EAR 

AREAS=AREAS+AREA2¥( XRATIOFSRATIC(CSI—-1LI#(SRATIO(C ISII/S20)) 
RCT( JJ) =2e* (AREAL +AREA2Z+ARFAS) 

DELTA( JJ) =ARFAI¥#FXRATIO%2 0/30 tAREA2K(DIST+eS*(XDIST#HDIST ) ) #tAREAS 
ALCAD(JJ)=(BMTOT (JJ) *#22*12¢)/DIST 
IF(EMTOT( JJ) -GTeXMPS)GC TC 779 

IFCABS(BMTOT(IJS)—-XMPS) eLE2e05)GC TO 780 

JJ=JI+1 

KJ=KJ+1 

GE 107.76 

CCATIAVE 

PFI=PKI-PINCR 

PINCR=e1*PINCR 

GC TO 776 

CCNTINUE 

CALL CALC 

Gc TO 11 

CCNTINUE 

CALL PLOT(0202020+9S9S) 

STOP 

EAD 


SLBROUTINE CALCI 


CCMMCKR ROT( 4000) 
COMMON JPLOT»LFAGE e TITLE (20) EL PLOT+SECT(10) »XMPSsADIST(4000) 


COMMCN AC 19) 0T(10) W610) 0E(10)sESK(10)+S(10)+SSH( 10) »«BMTOT( 4000) 
CCOMMCK BPHI (4000) eJJeK Js ALCAD (4000) sDELTA(4090) sAIX 

H=T(3) 

TT=T(2) 

DO=T(1)726 
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YEAR= (ACL ¥CHF2 e FTT40D) HAC 2) FC HEL SC S*TT) FAC 3) FC TT He SEH) HA( 4) * 05*TT) 
IS(CA(I)+AC 2) 4A 2)4A(4)) 
WRITE(Es100) YEAR 
100 FORMAT (1eXs*CENTROIC GCF SECTION IS *.,F1I0-3) 
AT X=A C2) ¥( TT ** 2/604 (HFtTTHTT/2 0-YEAR) €#*24+(TT/20—-YBAR) ¥¥2) tA(3) ¥(H¥EE 
Lesleot (TT +H/26-YEAR )*¥*¥2)4A(1)*( (0 2e¢*TT+HtEDD—-YSBAR) €*€2) 
RETURN 
ENO 


= 


SLBROLULTINE CALC 
CCMMCN fOT( 4000) 
COMMON JPLOT.LPAGEsTITLE(2C) -LPLUT.»SECT(10) »XMPS,ADIST( 4000) 
COMMCN A(19)+6T(10)5Ww(10)5E(10 )sESH(10)+S(10)+SSH( 10) +BMTOT( 4000) 
COMMON BFPHI (4000) sJJsKJsALCAD( 4000) »-DELTA(4000) »sAIX 
WEITE (6920) 
20 FCRNAT(12X¥,*MCMFRT CURVATURE RELATICNSHIPS -—- COMPOSITE fEAMS —? 
2* NEGATIVE BENDING® //) 
WRITE (6025) : 
2S FCRMAT ( U BEAM DATA (REBAR» FLANGE eWEBsFLANGEs AND COVERPL 


2 COVEFPLATE RESPECTIVELY )*// g AREA THICKNESS wID 
Sie e ESH SIGMA Y SIGMA SH*%9/) 
L=JJ-10 


CERO OVK= eS 
WRITE (6024) ACK) o TICK) oH (K)cECK) eESK(K) 9 S(K) » SSH(K) 
24 FCRMAT( 1X2 3F106¢3e2F1l2eir2Fl1l0e6 ) 
100 CCNTINUE 
WRITE (606246) AIX 
6246 FORMAT(12Xe"*MOMENT OF INERTIA OF SECTION 1S *®sF&e30*% INCHES*##4%5/) 
WRITE (626245) 
€245 FCRMAT( *CHECK® »//) 
DO 20C I[=L2JJ 
WEITE (69300) B8MTCT(1)s68PHI(I)seADIST( 1) sALOAD(T) -DELTACI) -ROT(I) 
BOO FCRMAT(1I2X% F100 eS sFl O0e6eFlOe3esF Be 3F10 06 92X oF l0 26) 
200 CCNTINUE 
CALL GRAPH 
RETURN 
ENO 


SLEROLTINE GRAPH 

CCMMCN FOT(4000) 

CCMMON JPLOT.LFAGE sTITLE (20) eLPLOTsSECT(10)5XMPS»*ADIST( 4000) 
COMMON A(10)eT(10)9WC1C) 9E(010)sESH(10) 5S (10) oSSH(10) sBMTOT( 4000) 
CCMMCN BPHI(4000) eJJx KU» ALGAD( 4000) sDELTA( 4000) sAIX 

TSCA Die talehetO COM SANO «JEL Gite Gillet) GAME EL ON Gl 2c — 1 9le19 — 3) 

IFCAC1) LT. 20001 DYPLOT=JPLOT+lL 

NNKK=0 
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BFHI(KJ+1)=0-0 
ALCAD(JSI+1)=0860 
RCTCIJ+1)=O064 
DELTA(JIJ+1)=064 
ACIST(JJt+1)=O06 
BMTOT(JI+1)=020 
TF C(LPLOT «EQ23)G0 
ITF (LPLCT2EQe4)G0 
IF (LPLOTeEQeS)GO 
IF(LFLOT.EQ26)G0 
IF(LPLOTeEQ.7)GO 
IF(LPLOT2EQe8)CO 
IF(LFLOT-EQ-9)GO 


IF (LPLOT-«EQe10)GO TO 40 
IF (LFLCTeEQe11)GO0 TO 41 
IF(LFLOTe-EQe12)GO TC 42 


BEHI{(KJ+2)=.004 
BMTOT(JJ+2)=606 
DELTA(JJ+2)=1e 
ALOAD(JJ+2)=306 
ACIST(JSJ#t2)=1 0 | 
PCT(IN+2)=203 
Ge Toss 
CCNTIAKUE 
DELTA(JJ+2)=le 
ALCAD(JJS#2)=156 
ADIST(CJJ+2)=1le 
RCT( JN +2)=203 
EMTGT(IS4+2)=2560 
BFEHI(KJ+2)=.0001 
GC TO 34 
CCNTINUE 
BNTOT(IIS+2) =S500 
PFHI(KI+2)=2.004 
DELTA(JSJ+2)=1le 
ALCAD(JJ+2)=306 
ARIST (J3t20=1 6 
ROT(JJ4+2)=063 
GC TO 34 
CCNTINUE 
BMTOT(JJS4+2)=2S56 
BEHI(KJ+2)=.0007 
DELTA(JJt2)=le 
ALOAD(JJS+2)=15S 0 
ACDIST(JSI#+2)=le 
RCT (JJ#2)=203 
GC TO 34 
CCNTINUE 

BMTOT (JJ#2)=256 
PFHI (KJ+2)=e0003 
DELTA(JJ+2)=025 
ALOAD(JJ+2)=15S 6 
ACES Utic)) = tte 
ROT( JS +2)=2005 
GC TO 34 
CONTINUE 
BPHI(KJ+2)=2002 
PNTCT(JSJ4+2)=6006 
DEL TACIU+2)=1 6 
ALOAC(JJ+2)=306 
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TO 
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TO 
TO 
TO 
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34 


ACIST(JJ+2)=1e 
ROT(JJ+2)=203 
GO TO 34 
CCNTINUE 
EMTOT(JJI+2)=606 
BPRI(KJ+2)=.0015 
DELTACIS#+2)=1e 
ALOCAD(JJS+2)=306 
ACIST(JJt2)=1 © 
ROT(JJ+2)=2-03 
GC TO 34 
CCNTINUE 
BMTCT(JJ+2)=606 
PPHI(KJ+2)=.-001 
DELTA(JJ#+2)=1le 
ALCAD(JJ+2)=300 
ACIST(JJ+2)=1le 
RCT(JJ+2)=.03. 
GO TG 34 
CCNTINUE , 
EMTCT(JJ+2)=S06 
EPHI(KJ+2)=.002 
DELTA (JJI#2)=16 
ALCAD(JJ#+2)=300 
ACIST(JJ+2)=1e 
RCT (JJ+2)=-03 
GC TO 34 
CCNTINUE 
BMTICT (JJ+2)=S06 
PPHI(KJ+2)=-001 
DELTA(JJ#2)=1.6 
ALCAD(JJS#2)=30-6 
ADIST(JIJt2) =1le 
ROT(JJ+2)=.03 

GO TC 34 
CCNTINUE 
EMTOT(JJ#+2)=S00 
EFRI(K J+2)=20005 
DELTA(JJ+2)=1le 
ALOAD(JJ+2)=306 
ACIST(JJ+2)=1le 
RiGihGI tei s03 
CCNTINUE 
IF(CA(1)eGTe-0001 
IF CAC 1) e-GTee 0001 
IF CAC 1) 2GT2 29001 
TF (A121) -GTe - 0001 


eANDeARAKK 0&Q00)CO0 
eANDeNKKK eEQe1)GC 
eANDeNARKK «EQ 22)G0 
eANDeNNKK cE O03) GO 


CALL 
CALL 
CALL 
GAEL 
CALL 


FLOT(0¢0 20002) 
FLOT(0e0e8e5Ss2) 
FLOT(11-028eSs2?) 
PLOT (11.020 e0e2) 
FLOT(0.000e002) 


e/Mefe Venklen ibs Bisa! 
CALL FLOT(1-2557 
GAL EEO Coe 7.Siei(, 
CALL FLOT(9-75 01 
GALE FEOT (1.2501 


CALL SYMBOL (10059 e0e00e1 eTITLES920430) 


02503) 
04592) 
045592) 
02572) 
02502) 


TO 
TO 
TO 
TO 


CALL SYMBOL (820 o207s0e1SeSECT 200043) 
CALL FLOT (2092e0s-3) 
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101 
102 
103 
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IF (NNKKeEQe1)GO TO 100 
IF(NNKKeEQe3)GC TO 105 
CALL AXIS(060:902001) 4HMCNENRT—FT eKIFSe 14050059000 sGMTST(JI+1) -BMTOT( 
1JJ4+2) +2020) 
IF (NNKK2EQ22)60 TO 104 
CALL AXIS(9.0006 e22HCURVATURE-RADIANS/ INCH s —22467 0080 00s 8PHI(KI+t1)o 
1BPHI(KJ+2)32020) 
25 CCNTIKUE 
IF (AC1)¢GTee C001 )CALL FLOT(-1465-170 0-3) 
CALL LINECBPHI»EMTOTsJJo10090) 
PECACID cle -OOOUM)GALE PROT (—2e6 1'Sies—3) 
NAKK=]1 
GO TO 34 
100 CCNTIAUE 
CALL AXIS(0e+e00es17HAPPLIED LCAD—K1IPS 01725029000 sALCAD(JJ+1) 2 ALOAD( 
1JJ+2)42020) 
CALL AXIS(066200¢+e29HCENTRE LINE DEFLECTION—-INCHESs—2937ex0 eeDELTACS 
1J4+1)2PELTA( JI+2) 520-0) , 
101. CCNTINUE 
TiCAVGS) iGo O O10 10) GALE ING 1.0 el) Ze s— 3) 
CALL LINEC(DELTAs ALCAD,JJe12020) 
EON DY OWS GOO EMER, MLR O24 SSG) 5 Ga) 
NAKK=2 
Gc To 34 
104 CCNTINUE 
CALL AXIS(0e200¢s LOGHROTATION—RADIANS 9-16 070 20 eeRCT(JJt1) eROT(IS 42) 5 
120¢) 
102 CCNTINUE 
IFCAC1) e-GTe e 0001 DCALL FLOT(14e2-1702—-3) 
CALL LCINECROTsBMTICTeJJ01 9000) 
IF CAC 1) eh Te 00001 )CALL FLOT(-20e01502-3) 
NNKK=2 
GC TO 34 
105 CCNTINUE 
CALL AXIS(OeeGee 1 4HMOMENT—FT eKIFS 2-14 050 200 eBMICT( JI+1 ) sBMTOT(JIt2 


1).,20.) 
CALL AXIS(Oes0ee28HNEUTRAL AXIS LCOCATION-INCHES s28e50290eeADIST (IS 
14+1).,ACIST(JJ+2).200) 
103 CONTINUE 
IF(AC1) eGTe eCCOl CALL FLOT (06. 0170-3) 
IF (A( 1) eGTee0001 CALL LINECBATOTsADISTeJJ5159020) 
0C 14 T=1+4000 
BMTOT(1I)=00 
PFHI(1)=06 
ALCAD(1)=064 
CELTA(T)=O6 
FCT(1)=00 
ADIST(I)=O6 
14 CCNTINUE 
FETURN 
END 
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